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Abstract. We consider a symmetric Anderson impurity model with a soft-gap 
hybridization vanishing at the Fermi level, Aj cx |a;|'" with r > 0. Three facets of 
the problem are examined. First the non- interacting limit, which despite its simplicity 
contains much physics relevant to the U > case: it exhibits both strong coupling 
(SC) states (for r < 1) and local moment states (for r > 1), with characteristic 
signatures in both spectral properties and thermodynamic functions. Second, we 
establish general conditions upon the interaction self-energy for the occurence of a 
SC state for U > 0. This leads to a pinning theorem, whereby the modified spectral 
function A{lu) — \uj\'^D{uj) is pinned at the Fermi level ut — for any U where a SC state 
obtains; it generalizes to arbitrary r the pinning condition upon D(lu — 0) familiar in 
the normal r = Anderson model. Finally, we consider explicitly spectral functions 
at the simplest level: second order perturbation theory in U, which we conclude is 
applicable for r < ^ and r > 1 but not for i < r < 1. Characteristic spectral features 
observed in numerical renormalization group calculations are thereby recovered, for 
both SC and LM phases; and for the SC state the modified spectral functions are 
found to contain a generalized Abrikosov-Suhl resonance exhibiting a characteristic 
low-energy Kondo scale with increasing interaction strength. 
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1. Introduction 

The Kondo effect, whereby an impurity spin is quenched by couphng to the low-energy 
excitations of a non-interacting metaUic host, has long occupied a central role in the 
study of magnetic impurities (see e.g. [1]). The effect is of course normally regarded 
as being dependent upon a metallic host, with low-energy and hence low-temperature 
impurity properties controlled by the non- vanishing host density of states at the Fermi 
level, a; = 0, and essentially independent of the details of host band structure. 

But what if the host exhibits semi-metallic character, with a spectrum whose low- 
energy behaviour exhibits a soft-gap at the Fermi level, phost oc with r > 0? 
There are quite a number of experimental candidates for such behaviour, ranging 
from semiconductors whose valence and conduction bands touch at the Fermi level 
[2], through heavy Fermion superconductors [3], to various two-dimensional systems 
including graphite sheets [4] and quasi-one-dimensional metals described by a Luttinger 
model [5]. 

The question above was first posed by Withoff and Fradkin [6], who studied the 

soft-gap Kondo model using a combination of 'poor man's' scaling and a large- mean- 
field theory (with N the impurity degeneracy). Since then, there has been much study 
of both the Kondo and the corresponding Anderson impurity models, in particular via 
scaling [6,7], large- A" expansions [6,8,9] and the numerical renormalization group (NRG) 
[7,10-13]. All these techniques, whether for the soft-gap Anderson or Kondo models, 
show the existence of two distinct types of ground state, between which in general a 
nontrivial zero-temperature phase transition occurs at a finite value of the host- impurity 
coupling (or, equivalently in the Anderson model, at a finite value of the impurity on- 
site interaction, U): a weak coupling or local moment (LM) state in which the impurity 
spin remains unquenched; and a strong coupling (SC) state in which a Kondo effect 
is manifest, and whose properties — in particular for the so-called symmetric strong 
coupling state considered here — have been argued to represent a natural generalization 
of Fermi liquid physics (see especially [13]). 

NRG studies in particular have devoted considerable attention to the spin-| (A^ = 2) 
particle- hole symmetric case, including thermodynamic [10-13] properties and, for the 
Anderson model, impurity spectral functions [12]. It is the symmetric spin-| soft-gap 
Anderson model we consider here, with aims that are modest, and threefold. First, 
to consider briefly the non-interacting limit, U — 0; second to establish rather general 
conditions upon the existence of a SC state at finite-^7; and finally, to examine the 
problem explicitly at the simplest possible level — straight second-order perturbation 
theory in U. 

There are two reasons for considering the non- interacting limit, simple though it is. 
First, and in contrast to the 'normal' Anderson model (r = 0), its behaviour is nontrivial 
and exemplifies much physics relevant to the interacting problem. An understanding of 
the impurity single-particle spectrum, Do{u!), is sufficient to study the non- interacting 
limit, since both 'excess' thermodynamic properties induced by addition of the impurity. 
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and local properties such as the impurity susceptibility, follow from a knowledge of it. As 
shown in section 3, both LM and SC states arise in the f/ = limit, with characteristic 
and distinct signatures in both spectral and thermodynamic functions, as indeed known 
in part from the work of ref. [13]. LM states alone are found to occur for r > 1 and 
SC states for r < 1, whereas finite-?/ NRG results [12,13] show that LM states alone 
occur for r > \ — an important contrast whose implications are considered in section 5. 
The second reason for considering the non-interacting limit is prosaic: a knowledge of 
it underpins finite-order perturbation theory in [/, as considered in section 5. 

In section 4, focusing on the finite-C/ impurity spectral function D{u!), we establish 
general conditions upon the interaction self-energy for the occurence of a SC state, 
from which follow in turn two results. First, that the low-frequency behaviour of the 
single-particle spectrum is D{uj) ~ 1^1*^, which is precisely the spectral signature of 
the SC state found in NRG studies [12]. Second, and relatedly, that interactions have 
no influence in renormalizing the low-tu asymptotic behaviour of D{lj). In consequence, 
one obtains a conservation on Ai^uj) = \u!\^D{u!) at the Fermi level, a; = 0: for any r 
where a SC state exists, A{u! — 0) is pinned at its non-interacting value for all U, a 
result that generalizes to arbitrary r the corresponding condition familiar for the r — 
Anderson model (see e.g.[l]). 

In contrast to the r = Anderson model where the predictions of perturbation 
theory in U about the non-interacting limit are well known (see e.g. [1]), the 
implications — and, indeed, general applicability — of a low-order perturbative treatment 
are not obvious for the soft-gap problem, and are considered in section 5 where we focus 
on the impurity spectral function. For r > 1 second-order perturbation theory is found 
to recover the characteristic low-o; spectral signature of the LM state found in NRG 
studies [12], viz D{u!) ~ |a;|^; and the resultant single-particle spectra are investigated 
in some detail. For r < | we find the SC state is indeed perturbatively stable upon 
increasing U from zero, and that the general conditions of section 4 for a SC state are 
satisfied at the second-order level; the resultant modified spectral functions A{u}) are 
also shown to bear a striking resemblance to that for the normal Anderson model, r = 0, 
exhibiting in particular the emergence with increasing f/ of a characteristic low-energy 
Kondo scale. For | < r < 1 by contrast, we argue that finite-order perturbation theroy 
in U about the non- interacting limit is simply inapplicable. 

There is a second, 'hidden' reason why we consider low-order perturbation theory: 
to illustrate its limitations, despite its strengths, even for r < | where the SC state 
is perturbatively continuable from the non-interacting limit. It is our belief that to 
describe analytically much of the underlying physics of the soft-gap Anderson model — 
and in particular to capture the transition between LM and SC (or generalized Fermi 
liquid) phases which render the problem of generic interest — requires, or at least invites, 
a new and inherently non-perturbative theoretical approach. We will turn to one such 
theory in a subsequent paper [14], for which the present work is in part a forerunner. 
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With the Fermi level taken as the origin of energy, the Hamiltonian for the spin-| 
Anderson model is given in standard notation by 

H H\iQgi -\- -ff impurity ~l~ -^hybridization (2.1a) 
= ^kflka + X](Q + ^fli^a)flia + ^ifel^Lcfea + cL^ja) (2.1b) 

fc,(7 fc,(T 

with Cfc the host dispersion, Vik the hybridization and the impurity level; for the 
symmetric case considered here, = — ^ with U the on-site interaction. 

We consider the zero-temperature single-particle impurity Green function, G{t) = 
—i{T{cia-{t)cl^}), with G{uj) expressible as 

G{uj) = [uj + irisgn{uj) - A(cu) - S(c<;)]"^ (2.2) 

where the limit rj — > 0-|- is henceforth understood. Here A(a;) is the hybridization 
function given by 

A{u;) = V '^^'r'' , , (2.3a) 

- efc + i?7sgn(cj) 

= Ar(cu) - isgn(cc^)Ai(cc;) (2.3b) 
and we consider throughout a symmetric hybridization 

A{u;) = -A(-a;) (2.4) 

(whose particular form is specified in §2.1). Prom particle-hole symmetry the Fermi level 

remains fixed at a; = V [/ > 0, whence the impurity charge n.j = ^^{nia) = 1\/ U; and 
the interaction self-energy S(a;) is defined to exclude the trivial Hartree contribution of 
{U/2)ni = which cancels = —U/2. S(a;) may likewise be decomposed as 

E(u;) = ER(a;) - isgn(u;)Ei(u;) (2.5) 

and the real parts of E or A follow directly from the Hilbert transform 

Fr(u;) = r ^Fi(a;OP i^A (2-6) 

J-oo TT \UJ — UJ\/ 

with F = S or A as appropriate. 

The single-particle impurity spectrum D{ijj) — — ;isgn(cc;)ImG'(ci;) follows directly 



as 



where 



and 



D{uj) = DP{uj) + D'^iuj) (2.7) 



D\.) = [A:(^) + E:(.)].-^ 

^ ' [u;-AR(a;)-ER(a;)f + [7^ + Ai(a;) + Ei(a;)]2 ^ ' 



D^iu) = 2. (2.8b) 

^ ' [^-AR(a;)-ER(a;)f + [77 + Ai(w) + Ei(^)]' ^ ' 
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Here D^{uj) refers to continuum (or 'band') excitations, while D'^{uj) allows for the 
possibility of discrete states reflected in pole contributions. 

The change in the density of states of the system due to addition of the impurity, 
Ap(c<j), will also prove central to the subsequent analysis. It is given by 



/\p{uj) = D{uj) 



UeGi.f-^ (2.9) 

TT OUJ 



did 

and is calculable directly once the impurity D{uj) is known. Note that equation 
while commonly derived explicitly for the non-interacting case (see eg [1]), is readily 
shown to hold generally for all U. From equation ( p.7|) , Ap{uj) likewise separates into 
'band' and 'pole' contributions, 

Ap(cu) = ApP(cj) + Ap^{uj) (2.10a) 

where in particular 



ApP(u;) = Z}P(cj) 



(2.10b) 



2.1. Hybridization function 

The hybridization function we consider, Ai{uj) (= Ai(— cu)), is given by 

= I ^« (M - r ■ i<\-\<oH (2.11) 

I : otherwise 

with r > 0. While a pure power-law hybridization is naturally not realistic on arbitrary 
energy scales, it captures in the simplest fashion the requisite low-a; behaviour; moreover, 
as for the r = Anderson model, one expects impurity properties to be controlled 
primarily by the low-u behaviour and to be largely independent of the details of host 
band structure. In general the hybridization thus contains a gap of magnitude 6, in 
which lies the Fermi level u; = 0; the 'normal' fiat-band Anderson model is recovered as 

a special case of equation ( p.ll| ) with r = = 6. There are four energy scales in the 

1 

problem, namely S, Aq"'' , D (the bandwidth) and U ; we choose to rescale in terms of 
1 

Aq""" , defining for later purposes 

cD = ^, D = -^, U = -^. (2.12) 

Ao^-'- Ao^-" Ao^-'- 

From the Hilbert transform equation ( p.6| ), AR(a;) = — Ar(— cj) is given by 

AR(a;) = duj, , ' .2 (2.13) 

where a principal value is henceforth understood. Notice from this that the wide-band 
limit D ^ oo, as commonly employed for the normal Anderson model r = = 6, can 
be taken only for r < 1. 
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For the gapped case, 5 > 0, we shall need solely the behaviour of AR(a;) for 
frequencies |u;| ^ | inside the gap. This is given from equation ( p.l3| ) by 



2An 



-uj- 



TT 



r-1 



B{5- r) + 



where B[5\ r) > is given by 

r2D/6 

B{6;r) = I dz- 



,2 • 



(2.14a) 



(2.14b) 



Our primary focus in sections 4 and 5 will be the 'soft gap' case, 6 = 0, where 
Ai{u) = Aq\u!\^. For this case, rescaling of equation (2.13) yields 



2 fD/l^ 



TT JO 



Ab,{uj) = sgn(u;)Ao|t^|''" 

which thus obeys the differential equation 
9Ar(o.) 



dy 



-Ar(u;) H 

UO TT 



d) 



(2.15) 



(2.16) 



Direct evaluation of equation ( [2.15| ) yields, for Ic^l < D, 



Ar(^) = -sgn(^)Ao jtan (f ) \uj\'- + 



vrr 



F(l,-r;l -r; \iu\/D) 



F{l,-r;l-r;-\uj\/D) 
2D 



-sgn(cu)Ao <^ tan (f ) 



\UJ\ 



7r(r -1) D 



+ 



\UJ\ 

u 



(2.17a) 
(2.17b) 



with F{a, z) a Gauss hypergeometric function. Equation (2.17) encompasses the 
results given in [13] and applies for any r > 0, including r = 1 for which the limit as 
r — > 1 of equation ( p.l7b|) gives 

3' 



AR(a;) = ^sgn(.;)|o;| In (M) + O 



\(jj\ 

u 



(2.18) 



showing the characteristic logarithmic behaviour that, as discussed in section 3, is 
indicative of the marginal nature of r = 1 for U = 0. 



3. Non-interacting limit 

As mentioned in section 1, the non-interacting problem is surprisingly rich: it contains 
already much underlying physics relevant to the interacting case and, for the soft-gap 
case in particular, gives rise to both SC and LM states as now considered. 

The impurity spectrum Do{uj) (with '0' referring to U = 0) is the primary quantity, 
since the excess density of states Apo(i^) follows from it via equation ( p.9| ). The low-u 
behaviour of the latter determines in turn the low temperature behaviour of the change 
in thermodynamic properties due to addition of the impurity; the 'excess' total uniform 
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spin susceptibility, specific heat and entropy being given trivially for non-interacting 
electrons by 



(T) 



2T .-oo 
2 /■°° 

T^J-oo 



du; Apo{u)f{u;) (1 - f{u)) 
dcu u;^Apo{uj)f{uj) (1 - fiiu' j) 



(3.1a) 
(3.1b) 



/oo 
dcu Apoicu) [fico) Inficu) + (1 - /(a;)) In (1 - /(a;))] (3.1c) 
-oo 

where f{uj) is the Fermi function (with chemical potential /i = for all T due to 
particle-hole symmetry), and ks = 1 has been taken. 

The band/continuum part of Apo(^^) is given from equations (2.9,10) by 



1 - 



du 



(3.2) 



with Dq{uj) from Eq(2.8a). For the gapless case 5 = 0, a straightforward calculation 
using Eq(2.16) for dA-^/dui gives a simple relation between ApQ(a;) and Dq{uj) for 
|a;| < D: 



:i-r) 

g(a;) 







where 



2An/^''"^ 



(3.3a) 



(3.3b) 



The pole contributions, DI{lo) and Apo(<^), are considered below. Here we simply note 
that for finite bandwidth, D, there are always such contributions outside the band, 
\uj\ > D. These however are of no importance to the problem, and are not considered 
explicitly in what follows where \u!\ < D is implicit. 



3.1. Gapped case: LM state 

We consider first the case of an insulating host, with S > 0. Since Ai(a;) 
< S/2 inside the gap, it follows from equation (2.8b) (with E = 0) that 

L>P(u;) = qSiu) 

with poleweight q given by 

fdA^{u)\ 

' w=0- 



1 - 



\ du 



1 + 



2Ao 



TT 



r-1 



B{S;r) 



for 



(3.4) 



(3.5a) 



(3.5b) 



where equation (2.14) is used. Prom equation (2.8b) the band contribution Dq{u) 
naturally vanishes inside the gap, and behaves as Dq{u) oc Ai(a;) ~ (|a;| —S/2Y close to 
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the gap edges. Dq{uj) is thus dominated by the discrete state at the Fermi level, a; = 0, 
as too is Apo(<^): from equations (2.10b) and (3.4,5), 

ApP = 5{uj) (3.6) 

whose poleweight of unity reflects a 'whole' single extra state at the Fermi level induced 
by addition of the impurity. 

These features — a whole excess state at the Fermi level, with non- vanishing weight 
on the impurity — are the hallmark of the LM state for [/ = 0. They naturally control 
the low-T excess thermodynamic functions, given from equations (3.1) by 

8kT 



XL,{T) = = ixcu.e(T) (3.7a) 



5Lp(0)=ln4 (3.7b) 

with corrections that are thermally activated (as too is Cj'^p(T)). The residual entropy 
of In 4 is physically obvious: the excess Fermi level state has four occupancies — empty, 
t-spin or J,-spin occupied and doubly occupied — with equal a priori probabilities for 
U = 0. This is also why Xi'inp(T') is half a Curie law: only half of the four occupancies, 
the singly occupied states, are paramagnctically active. This situation will of course 
change immediately for any U > where doubly occupied (and hence empty) spin 
configurations are suppressed, producing instead XimpiT) = Xcurie(2^) as T — >• 0, and 
^imp(O) = ln2. 

The gapless case, (5 = 0, is considered in sections 3.2,3, but first we ask what 
happens as 5 — > 0. From equation (2.14), 

B{S; r) = -— - + -— - F(l, -r; 2 - r; -x) (3.8a) 

sm(7rr) [r — 1) [1 + x) 

with 

x=^. (3.8b) 

From equation (3.5b), the behaviour of the poleweight q for 5/2D <S 1 can thus be 
obtained. For r > 1 one finds 

7r(r — 1) 

which remains finite as 5 ^ (where the resultant q = q{uj = 0), from equation (3.3b)), 
showing the persistence of the LM state for r > 1 when 5 = 0. For r < 1 by contrast, 

+ 0(1) : r<l (3.10) 



1 + 



+ 0{S'-^;S) :r>l (3.9) 



2r <- .-1 

= 1 + -^AoD^ 
sm 7rr 



2D_ 

which diverges as 5 — (as too does the marginal case of r = 1, where 

g"^ ~ l/ln{2D/S)). For r < 1, the LM state does not therefore survive closure of the 

gap. 
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(a) 




(0 = 



(0 = 



(0 = 



Figure 1. Schematic of non-interacting impurity spectra Do{ui) for (a) LM state, 
r > 2; (b) LM state, 1< r < 2; (c) SC state, r < 1. 



3.2. Gapless case: r > 1, LM state 

For r > 1, the LM state indeed persists with the gap shut, 5 — since Ai(a;) ~ 
decays to zero as cj ^ more rapidly than Ar(cj) ~ 0(|cj|) (see equation (2.17b)), 
equation (2.8b) yields directly = qS{uj) with q given by equation (3.5a) and hence 
equation (3.9) with 6 = 0. Likewise, from equation (2.10b), Apo(a;) = S{u!) contains a 
'whole' extra state. 

In contrast to the gapped case however, Dq{lj) now extends down to a; = 0. From 
equations (2.8a) and (2.17b) one obtains 



\U) 



r-2 



TT 



+ o(kl ) 



(3.11) 



diverging as a; — > for 1 < r < 2 and vanishing for r > 2, as found in [13] . The full Do{u;) 
consists of course of both pole and band contributions; it is illustrated schematically in 
figure 1. From equations (3.3) and (3.11) the low-o; behaviour of Apo(a;) is 

Aog, 



r)\uj\ 



ir-2 



+ 



TT 



(3.12) 



The excess thermodynamic properties follow directly via equation (3.1): 'S'i'^p(O) = ln4 
again, while 

xLp(r) = iXcuHe(T) + c'T^-' + 0{T'^-'; (3.13a) 



(T) = d'r-^ + ©(T^^'-^); T^+') 



(3.13b) 



where c' and d' arc negative constants (whose sign reflects the low-cj 'depletion' of 
Apq{uj) — the coefficient of the leading |c<j|''^^ term in equation (3.12) being negative 
for r > 1). Not surprisingly, and in part for the reasons discussed in section 3.1, the 
behaviour equation (3.13) characteristic of the U = limit differs from that found by 
NRG for the [/ > LM regime [12]: Ximp(T) = XcuHe(T) + c'^T^'^ and amp(T) oc T\ 
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3.3. Gapless case: r < 1, SC state 



Prom equations (2.8a) and (2.17b), a simple calculation gives the low-o; behaviour of 
the continuum contribution to Do{u!), viz 

/? 2 



\UJ\ 



il-r 



+ 0(|a;|'-''') (3.14a) 



I + /32 Aog(O) 
where 

/3 = tan(fr). (3.14b) 
Apq{uj) follows directly from equation (3.3a) and has the same leading asymptotics as 



Now consider Dq{uj), given from equation (2.8b) with S = 0. From equation (2.17b) 
the leading low-o; behaviour of Aji{iu) for < r < 1 may be cast as 

AR(a;) = -sgn(a;) tan (f r) Ai(a;) 

whence 

L>o» = L(Ai;r) 



(3.15) 
(3.16) 



with L{x; A) (used again in section 4) defined generally by 

L{x]X) 



rjTT 



tan (f a) 
A 



X 



+ [r] + x] 



77^0 + 



6{x). 



Since Ai(a;) 



tan (f a) 
Ao|a;|'', it follows that 



1 



\UJ 



1-r 



5{u;) 



= 



(3.17a) 
(3.17b) 

(3.18) 



Ao tan (f r) 

i.e. there is no pole contribution to Dq{uj) itself, consistent with section (3.1) as 5 0. 
Hence Do{u) = D};{ u) is given by equation (3.14), with Dq[u!) ~ |a;| ^ as found in [13] 
(and illustrated schematically in figure Ic). 

But from equation (2.10b), Apq{uj) is given in contrast by 

ApKu) = D^{uj) [1 + Aortan (f r) jc^r"^] (3.19a) 
= rS{u;) (3.19b) 

(via equation (3.18)), and thus contains a 5-function contribution of weight r, as well 
known from the work of [10,13] where the result was obtained from analysis of the phase 
shift. This behaviour is the characteristic of the SC state, and suggests the interpretation 
[13] that a fraction r of a conduction electron occupies the decoupled excess state at the 
Fermi level, the remaining fraction 1 — r being absorbed into the resonant continuum 
centred on a; = 0. And note again that, in contrast to the LM state, the excess level 
has no weight on the impurity itself: D^{uj) = 0. 
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Excess thermodynamic properties follow directly via equation (3.1), namely 

xLp(r) = iXcuHe(T) + cT- + 0{T'-'^) (3.20a) 
C^^^{T) = dT'-' + 0(T2(i-'')) (3.20b) 
'^imp(^) = 2r In 2 + eT^'"^ + 0(7^^^-'')) (3.20c) 

where c, d and e are positive constants. These are worth noting for comparison to NRG 
results obtained for the C/ > SC phase, namely [10-13] 

Ximp(T) = iXc.He(T) + c[T-' + c'^T-^' (3.21a) 

Ci„,p(T) oc T'-' (3.21b) 

5i„p(T) = 2r In 2 + e'T^'^ (3.21c) 

It is well known [10,13] that the leading low-T behaviour of ximp and S'lmp for the SC 
phase — viz |xcurie and 2rln2 respectively — are given precisely by the U = limit result. 
But it is striking to note that the leading T-dependences of Cimp(T) and AS'imp(T), 
namely T^~'^, are also inherent to the non- interacting limit, as too is the T"^ correction 
to Ximp- Only the T contribution to the NRG Ximp — ^which applies only for ^ < | 
where the SC phase arises for U > — is absent in the non-interacting limit. 

3.4- Local impurity susceptibility 

The 'excess' Xi'mp(^) discussed above refers to the change in the total uniform spin 
susceptibility of the system induced by addition of the impurity. Here we comment 
briefly on the local impurity susceptibility XuiT), which has also been studied via NRG 
[10,13]. It is defined by 



Xii{T) = -g^B- 



dh 



(3.22a) 

h=0 



with Siz referring to the impurity spin and h a magnetic field acting solely on the 
impurity; and is given in standard notation by 

Xii{T) = {g^B? f dr(^,,(r)^,,) (3.22b) 

with (3 = 1/T. In the non-interacting limit, x'iiiT) is trivially evaluated and given by 

Z J—oo J —CO U}\ — 102 

Separating Dq{u}) = q5{uj) + D^iuj), and using particle-hole symmetry, a simple 
calculation gives 

{giiBf ST 2i-oo oj V2r 

+ i d.. r DlMDt(..)ll^^2l^lM\, (3.24) 

Z J—oo J—oo OJi — OJ2 
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For the gapless LM regime r > 1, where the poleweight q{r) 7^ is given by equation 
(3.9) (with (5 = 0), the low-T behaviour of xl{T) is thus 

Xlm = ^XcuHe(T) : r > 1 (3.25) 

with leading corrections O [min(T''~^, 1)] arising from the second term of equation (3.24). 
Hence, as expected for a LM state and mirrored also in x°iiip(T') equation (3.13a), the 
impurity spin remains unquenched; although X°(T')/x?jjp(r) = g^<lasT'— >0, 
reflecting the fact that the 'whole' excess level induced by the impurity has only partial 
weight on the impurity itself. 

For the SC phase r < 1, by contrast, the excess level has no overlap on the impurity, 
g = 0. Only the final term in equation (3.24) survives, and using equation (3.14a) is 
0(1) for r < ^ and 0(T^~^^) for ^ < r < 1; hence, in particular, 

limrx°(r) = : 0<r<l. (3.26) 

As for the normal (r = 0) Anderson model, the impurity spin is thus locally quenched 
in the entire SC regime r < 1; in contrast, for obvious reasons, to the behaviour of Xhnp; 
equation (3.20a). 

The above behaviour — complete (SC) versus incomplete (LM) quenching of the 
impurity spin — is also found in NRG studies of the soft-gap Kondo model [11,13]; and 
the total spin quenching symptomatic of the SC phase is one reason why it may be 
regarded [13] as a natural generalization of conventional Fermi liquid physics. 

We have seen that even the non-interacting limit contains both LM and SC states, 
whose characteristics are reflected in the behaviour of the impurity spectrum Dq{uj) — 
by whether (LM) or not (SC) there is a (5-function contribution at the Fermi level — and 
hence in turn in Apo(^^) and resultant thermodynamic functions. One final point should 
however be noted: the f/ = 'phase diagram' consists of LM states for r > 1 and SC 
states for r < 1. This is in contrast to what is found by NRG for f/ > [12,13], where 
SC states arise only for r < | and where only LM states occur for r > |. We consider 
the implications of this further in section 5. 

4. U > 0: conditions for SC phase 

For non-vanishing interaction strengths, f/, we now seek conditions under which, upon 
increasing U from zero, a SC state will persist; (only the gapless problem, and for r < 1, 
need be considered, since for r > 1 the U = Q ground state is a local moment one). 
Employing a phase shift analysis that parallels Nozieres' Fermi liquid description of the 
'normal' Kondo effect [15], Chen and Jayaprakash [10] have argued that the SC state is 
characterized generally by 

ApP(cu) = r5{oj) : SC (4.1) 

— precisely as for the non-interacting limit. 
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The question now is: under what conditions upon the self-energy S(a;) will this 
arise? Ap^iu) is given by equation (2.10b) with D^{uj) from equation (2.8b). Since Aj 
and Si are non-negative and, for |c<j| < ^lioo) = Aq\uj\^ vanishes only for u; = 0, 
the only 5-funtion contribution that could arise in ApP(a;) for \uj\ < D is of course for 
a; = 0. Hence only the low-a; behaviour of A and S is relevant, and ApP(a;) is given by 
(of equation(3.19a)) 

ApP(cu) = DP{uj) [l + Aortan (f r) |cu|"-^] (4.2a) 

with 

DP(u;) = 2 (4.2b) 

^ ^ [AR(cc;) + SR(c.)f + [^ + Ai(a;) + Ei(u;)]' ^ ' 

(where the 'bare' u contribution to equation (2.8b) can again be neglected since, for 
r < 1, it is subdominant to An{u!) ~ Notice also that if the low-a; behaviour of 

Ei(a;) is of form 

Ei(t^) = a\uj\^ : u^O (4.3a) 

with — 1 < A < 1, where a = a{U) is a (necessarily positive) constant, then the 
corresponding low-o; behaviour of ^^(a;) follows directly from the Hilbert transform 
equation (2.6) as 

Er(u;) = -sgn(a;) tan (f A) Ei(u;) : u ^ 0. (4.3b) 

If Ei(a;) and hence ER(a;) decay to zero as a; — > more rapidly than Ai(a;) — i.e. 
if Ei(a;) is of form equation (4.3a) with r < A — ^then, trivially, the low-o; behaviour of 
equation (4.2b) for D^{uj) is precisely that of the U = limit, i.e. Dp(cj) = L(Ai;r) 
(see equation (3.16)). Hence, as in equations (3.17) ff, ApP(a;) = rS{uj) arises. If by 
contrast Ej/r dominate the low-cu behaviour of equation (4.2b) — i.e. if Ei(c<j) is of form 
equation (4.3a) with A < r — then from equations (4.3) and (3.17), D^{uj) = L(Ei;A); 
from equation (4.2a) it then follows that ApP(a;) ~ |a;|''~'*'5(a;) = since A < r, i.e. there 
is no pole contribution to Ap{uj). Finally, if the low-a; behaviour of Ei/rIs the same as 
that of Ai/R, viz A = r in equations (4.3), then a directly analogous calculation gives 
ApP(a;) = r/[l + a{U) / Aq] S{u;); i.e. a 5-function contribution but with a [/-dependent 
weight that is less than r. 

The question posed above is thus answered: for ApP(a') = rS{uj) to arise, and thus 
a SC state to be realized for U > 0, Ei(a;) and hence ER(a;) must decay to zero as a; — > 
more rapidly that |a;|'', i.e. 

Ei(a;) Q;|a;|^ : r < A (4.4) 

This has important implications for the low-a; behaviour of the impurity single- 
particle spectrum D{u;) (= D^{u;)), since from equation (2.8a) it follows that the low-a; 
asymptotic behaviour of D{u!) is precisely that of Do{u!), whence in particular 

D{uj) la;]-'- (4.5) 
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as indeed found for the SC phase in NRG studies of the impurity spectrum [12]. More 
importantly, using 

hm br^fo;) = hm U''Do(u;) (4.6) 

together with equation (3.14), and defining A{uj) = |c<;|'"Z^(c<;), it follows that 

ttAo [l + tan=^ (f r)] A{u; = 0) = 1. (4.7) 

We believe this to be significant. For the 'normal' Anderson model, r = 0, it recovers the 
well known result (see e.g.[l]), normally viewed as a consequence of the Friedel sum rule, 
that the impurity single particle spectrum is pinned at the Fermi level uj = 0, i.e. that 
D{u! = 0) = I/ttAq (in this case for all U). Equation (4.7) represents a generalization 
of this pinning condition to arbitrary r where a SC state obtains, whose continuity in 
r reflects the fact that interactions have no influence in renormalizing the asymptotic 
behaviour of D{u!) as a; — > 0; and which is entirely consistent with the conlusions of 
Gonzalez-Buxton and Ingersent [13] from NRG studies that the SC state embodies a 
natural generalization of standard Fermi liquid physics. The extent to which equation 
(4.7) is captured in practice should also provide a good test of the accuracy of NRG 
calculations at low frequencies, and will be discussed elsewhere; moreover the generalized 
pinning condition will prove central to our local moment approach to the problem, as 
will be discussed in a subsequent paper [14]. 

One important question is not of course answered by the above considerations: 
for what range of r will the condition equation (4.4) for a SC state actually arise? 
NRG calculations give r < | for the SC state [10-13]. We examine this question in 
the following section within the framework of second order PT in U, together with a 
corresponding analysis of the evolution with U of the single-particle spectra appropriate 
to the LM state. 



5. Perturbation theory in U 

Low order perturbation theory in U about the non-interacting limit is probably the 
simplest and certainly the most conventional approach to the problem. For the 
'normal' Anderson model, r — 0, its predictions are of course well known (see e.g. [1]): 
while restricted by contruction to weak coupling interactions U, and thus incapable of 
capturing strong coupling 'Kondo' asymptotics, it generates order by order characteristic 
Fermi liquid behaviour, in particular that Si ~ 0{uj'^) and that the impurity spectrum 
is pinned at the Fermi level, D{0) = I/ttAq; and the single-particle spectrum evolves 
continuously upon increasing U from the non-interacting limit, in accordance with the 
fact that the normal Anderson model is a Fermi hquid for all U/Aq > 0. For r > by 
contrast the implications of a low order perturbative treatment — and even whether such 
is in general applicable — are far from obvious, and are considered here at the simplest 
second order level. 
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Figure 2. Second order self-energy diagram, witii bare {U = 0) impurity propagators 
denoted by solid lines and the on-site impurity U by wavy lines. 



The familiar second order self-energy diagram is shown in figure 2 and may be 
written as 

E(a;) ^U' r ^ °n(a;i)Go(c^i + u;) (5.1) 

J-oo 2m 

where Gq{lj) is the non-interacting impurity Green function with spectral representation 
G„H = r d.. £pL^^, (5.2) 

The 'polarization bubble' ^{u) is given by 

°n(u;) = i r ^GoMGoioui - u) (5.3a) 

= (5.3b) 

(where the latter follows from a trivial change of variables in equation (5.3a)); its spectral 
representation is 

r°° dcui Im °Q(a;i)sgn(a;i) 
cui — cu — ir)Sgn{u!) 

Prom equation (5.1), using (5.2) and (5.4) together with Do{u;) — Do{—u;) (particle- 
hole symmetry), a straightforward calculation gives 

Eiioo) = U'^ /''^' dooilm ^{u;i)Do{u;i - \u;\) (5.5) 
Jo 

which is thus readily calculable from a knowledge of Dq{uj) (section 3) and ''n(a;). 
Similarly, using equation (5.2) again, a directly analogous calculation yields 

^Im 'tl(a;) = /''^'da;iDo(^^i)^o(^^i - l^^l) > 0. (5.6) 
^ Jo 

Note that, using equation (5.6), equation (5.5) implies Si(a;) > as required by 
analyticity. Equations (5.5) and (5.6) are the basic equations to be analyzed, as now 
considered for the gapless problem (S — 0). 

Separating Dq{u!) into pole and band contributions, and using Dq{u!) — q5{u!), 
equation (5.6) reduces for < to 

^Im °n(a;) = lq^6{u) + qD];{u) + du;iZ^o'(^i)^o('^i - kl)- (5-7) 

t/ 



= r dc^ Im ^(a;Osgn(c.O 
J-oo TT UJt — UJ — insRuiuj) 
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This encompasses both r > 1 where g 7^ is given from equation (3.9) with 5 = 0; 
and r < 1 where g = (section 3.3) and only the 'band-band' contribution to Im ^{u) 
survives. [We add in passing that pole-contributions to Dq{uj) from outside the band, 



> D, generate spectral contributions to Im 'll(c<j) for D < \uj\ < 2D; as we shall be 
concerned only with the low-a; behaviour of Im ''n(a;) and hence Si(a;), we refrain from 
showing these, although they are fully included in numerical calculations where they are 
necessary to ensure the correct normalization of -D(a;).] 

The low-o; asymptotics of Im °Q(a;) follow from equation (5.7) using the results of 
section 3 for Dq{u;). For r > 1 we obtain 

^Im°n(o;)"~°|g2^(a;) + gD^(o;) + 0(|a;|2'-3) : r>l (5.8) 

where the 0(|a;p''~^) corrections arise from the band-band piece of equation (5.7) 
and diverge less rapidly as a; ^ than the leading \uj\^~'^ behaviour of Dq{u) 
(equation (3.11)). For r < 1 by contrast, only the band-band contribution to equation 
(5.7) survives, and a simple calculation using equation (3.14) gives the leading low-o; 
behaviour 

^Im'tl(a;)''~°C(r)|a;|^-''- : r<l (5.9a) 



with C > given by 



cos^ I f r ) /•! 



C{r) = ^ / dy y-^{l - (5.9b) 



cos^ (|r) 2^'""^ r(l 



Ao r( 



(5.9c) 



These results may now be used in equation (5.5) to determine the crucial low-o; 
asymptotics of Ei(a;), as now considered separately for r > 1 and r < 1. 

5.1. r>l: LM state 

From equations (5.5) and (5.8) the low-o; behaviour of Ei(a;) is given by 

Ei(a;) = --^ \q5{u) + ^.D'^iu) + 0{\u\^'~'')\ (5.10a) 



4 



2^3 



^ 5{uj) + Y.\{uj) (5.10b) 

where, using equation (3.11), the asymptotic behaviour of the 'band' contribution Sj (0;) 
is 

The low-a; behaviour of Si(a;) is thus integrably singular, and manifestly not that of a 
Fermi liquid or any natural generalization thereof: not surprisingly, since the underlying 
[/ = ground state for r > 1 is a local moment one. The corresponding real part, ER(a;), 
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follows from the Hilbert transform equation (2.6) using equation (5.10a), and has the 
leading low-o; behaviour 

SK(..)^-°^pf-) (5.11) 



4 Ku, 

(with corrections 0(|u;|''~^, |a;|) arising from the transform of E^(u;)). 

From equation (2.8) the low-a; behaviour of the impurity spectrum D{u!) thus 
follows directly as 

This behaviour — D{uj) ~ \uj\^ — is precisely the spectral hallmark of the LM regime 
found in NRG calculations for U > [12]. It is of course in marked contrast to what 
obtains in the non-interacting limit, viz equation (3.11) -Dq(ci;) ~ I^^T"^ (together with 
a pole contribution for [/ = 0, which is ehminated entirely for U > 0). Although the 
lowest frequency spectral asymptotics thus change character abruptly upon increasing U 
from zero, it is however straightforward to show (using the asymptotic forms of Sj^/E^ 
above) that for sufficiently small U there exists a crossover scale uq = \Uq^ such that 
for [D S>)|a;| ^ ojq the behaviour of D{(jj) is that of the non-interacting limit, viz 
Diuj) ~ while for |u;| <^ ujq, D{uj) ~ \uj\^ as in equation (5.12). 

Representative single-particle spectra for the LM state, obtained at the second- 
order level, are illustrated in figures 3 and 4 for r = 1.5. We consider first a 'weak 

1 

hybridization' example which, for r > 1, entails D = Aq'^ D <^ 1 (see equation 

1 

(2.12)). For D = 1.5 X 10^^', figure 3a shows the dimensionless D'{ijj) = Aq' D{uj) 

vs. uj = uj/Aq'"^ for three reduced interaction strengths U = U/Aq~'' = 2.5 x 10~^, 
5 X 10^^ and 7.5 x 10~^. The dominant visible feature of the spectra are the Hubbard 
satellites, which for all C/'s shown are centred to high accuracy on cJq = ^Uq^ (with the 
poleweight q 0.95 from equation (3.9) with 5 — 0). This is expected physically: the 
weak hybridization regime is 'close' to the atomic limit, Aq = = Vik (where q — 1), 
which by a well known accident for the particle- hole symmetric case (see e.g. [1]) is 
captured exactly by second order PT, and where equation (5.11) with g = 1 is exact for 
all UJ. The sharp Hubbard satellites of figure 3a thus correspond simply to weak resonant 
broadening, and only slight shifting, of the atomic limit poles occuring at uj = ±|f/. 
The low-o; behaviour of the spectrum, equation (5.12), is not directly visible in figure 
3a, but is clear from the inset to figure 3a and in figure 3b. In the latter, for the same 
parameters as figure 3a, we show log {TT[Uq]'^ / 12) D'{uj) vs. logo). That D'{uj) ~ |a;|'' 
as a; —> is evident; as too is the accuracy of equation (5.12) in its entirety. The above 
mentioned crossover to |a;|^~^ behaviour for |a;| ^ (Do is also evident in figure 3b for the 
lowest U example. 

The spectral features naturally evolve smoothly with increasing hybridization 
strength. For r = 1.5 again, figure 4 shows D'{uj) vs. uj for a 'strong hybridization' 
example, D = 10 (where q ^ 0.2 <^ 1), for [/ = 4, 8 and 20. The low-cD behaviour, 
D'{uj) ~ |a;|'', is clearly seen in all cases. The Hubbard satellites are again centred on 
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Figure 3. LM regime, r = 1.5. (a)L)'(tD) = Aq~'' D{lo) versus (D — lu/Aq^'' for 

reduced interaction strengths U = 2.5 x 10""' ( ), 5 x 10""' ( ) and 7.5 x lO"** 

( ); all for reduced bandwidth .D = 1.5 x 10~^. Inset: low-(D behaviour of D'{lj). 

(b) iriUqf' /12 D'{u}) versus cD on a log-log scale for same parameters as (a). 



ljq = \Uq^ for sufficiently low U (e.g. for JJ = 4); and in all cases, since q <^ 1, lie well 
below the scale of |f/ characteristic of the atomic limit. It is also seen that the satellites 
become increasingly diffuse with increasing interaction strength, although whether this 
is a genuine feature is not clear since second order PT is by construction confined to 
weak coupling. 

5.2. < r < |.- SC state 

For r < 1, Im tl(a;) is given by equation (5.9) and the a; — > behaviour of 
Dq{uj){= Dq{uj)) by equation (3.14). Hence from equation (5.5) the low frequency 
behaviour of Ei(a;) is 

Si(a;)"~M(t/)|u;|2-3'^ (5.13a) 



Magnetic impurities in gapless Fermi systems: perturbation theory 



19 



4.0 



3^ 

Q 2.0 



0.0 



-1.0 



-0.5 




OX) 



0.5 



1.0 



Figure 4. LM regime, r = 1.5 and reduced bandwidth D = 10. D'{Cb) versus cD for U 
= 4 ( ), 8 ( ) and 20 ( ). 



with 



cos^ ( —r 1 /■! 
(U) = ^ C{r) / dy y'-'^{l - y)' 

Aq Jo 

^(H^(^)r(2[l-r])r(l-r) 



.COS 



Ao 



r(3[l-r]) 



(5.13b) 



and C(r) given by equation (5.9). From the Hilbert transform equation (2.6), the 
corresponding behaviour of XlR(a;) follows as 



SR(a;) 



-sgn(u;)7|a; 
— sgn(a;) tan 



(5.14a) 



1(2 -3r) a{U)\uj 



|2-3r- 



0<r<| 
r>\ (5.14b) 



(with logarithmic behaviour ~ |cc;| In |cc;| for r = |), where 7r7 = da; Si(a;)/a;^ > 0. 

The important point here is that for r < \, where r < 2 — 3r, Si(a;) and ER(a;) 
decay to zero as a; — > more rapidly than Aj/r ~ ■ The conditions established in 
section 4 for the SC state are thus satisfied at the level of second order perturbation 
theory, and the SC state is hence perturbatively continuable from the non-interacting 
limit. Interactions do not renormalize the asymptotic behaviour of D{uo) as a; — 0, and 
the generalized pinning condition equation (4.7) is thus satisfied. Note moreover that, 
as for the normal Anderson model r = 0, Si(a;) vanishes at the Fermi level a; = 0. 

Representative single-particle spectra for the SC phase are shown in figure 5 which, 
for r = I and the wide-band hmit {D — oo), shows D'{Cj) vs. o) for [/ = 1, 5 and 
7.5. As found in NRG studies [12], the gross features of the spectrum are dominated by 
the low-o; behaviour where D{ijj) ~ ■ And with increasing interaction strength [/, 
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the low frequency resonant continuum narrows and emergent Hubbard sateUites become 
evident as weak shoulders in the spectrum. 



1.0 



^ 0.5 




Figure 5. SC phase, r = \. D'{lo) = A^-'' D{lo) versus w for C/ = 1 ( ), 5 ( ) 

and 7.5 ( ). 



A far more revealing expose is seen in figure 6 where, for the same parameters as 
figure 5, ttAq 1 + tan^ (^|rj A{uj) vs. a) is shown, with 

A(u;) = \u;\''D(u;). (5.15) 

Note first that the generalized pinning condition equation (4.7), which applies for 
any r where a SC state obtains, is manifestly satisfied. The low-o; behaviour of 
A(u}) is cusp-like, and from the asymptotics given above it is readily shown that 



where the coefficient ai oc tan(7rr/2). For 



l-r 



while for 



< r < 2 it behaves as 



A{uj) ~ A{0) - [ai|cj|i-" + 021^^12(1-20 
< r < i the cusp is thus of form AA{uj) ~ \uj 
AA(uj) ~ |a;|2(i~2'-). g^j^j f^j. 7- = 0, the parabolic behaviour AA{u;) ~ cj^ characteristic 
of 'normal' Fermi liquid behaviour is recovered. In fact, the behaviour of A{lj) shown 
in figure 5 for r = 0.25 is strongly reminiscent of spectra characteristic of the normal 
Anderson model, r — 0; the latter being shown in figure 6 for the same U values. 
The parallels are obvious, in each case A{u; — 0) being [/-independent (pinned) and 
with Hubbard satellite peaks progressively evolving with increasing interaction strength. 
Most significantly, we see in either case the emergence with increasing [/ of a low-energy 
scale, reflected in the half-width, ujk, of A{uj) which narrows progressively as U is 
increased. For r = this is just the emergence of the ordinary Kondo scale, while for 
< r < ^ it is the generalization thereof known in a thermodynamic context for the 
soft-gap Kondo model itself [6,10]. 




We shall not however pursue here the evolution of the Kondo scale with increasing 
U since the above analysis, while showing that the SC state is pcrturbatively continuable 
from U = and enabling a description of the spectrum at low U, also points clearly to the 
limitations of a low-order perturbative treatment. For the normal Anderson model, we 
know that low-order perturbation theory is intrinsically incapable of describing strong 
coupling 'Kondo' asymptotics — D{lj) narrowing algebraically with increasing U and 
failing thereby to yield the correct exponentially small Kondo scale cuk ~ exp(— ttC/ /SAq) 
(see e.g. [1]). And the situation is even more acute for the SC phase with r > 0: here, 
from NRG studies of the gapless Anderson model [12,13], it is known that for < r < | 
there is a critical line Uc{r) below which a SC state obtains, and above which the 
ground state is by contast a LM one. As C/ — > Uc{r)—, the low energy Kondo scale 
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ujk characteristic of the SC state must vanish, and this clearly cannot be captured 
by straight perturbation theory. To capture such behaviour, and hence in particular 
to describe the SC/LM phase boundary, an intrinsically non-perturbative approach is 
needed, and will be given in a subsequent paper [14] . 

5.3. \<r <l 

Finally, we turn to the regime | < r < 1. Here the low-u; asymptotics of the second 
order self-energy are again precisely those given in section 5.2, viz equations (5.13) and 
(5.14b). But in contrast to < r < |, Si(a;) and SR(a;) do not decay to zero as o; — >• 
more rapidly than Ai/r ~ |a;|'', whence the conditions established in section 4 for the 
SC state are not satisfied for | < r < 1: there is no pole contribution to Ap{uj), and 
the pinning condition equation (4.7) is not satisfied since the low-o; behaviour of the 
impurity spectrum is apparently dominated by Sj/r, viz D{u!) ~ [Ei(a;)]~^ ~ — 
diverging as Ct; — > for ^ < r < |, but less rapidly than the behaviour characteristic 
of the SC state; and vanishing for | < r < 1, but less rapidly than the \uj\^ behaviour 
symptomatic of the LM state. Superficially, therefore, one appears confronted by a state 
that is neither SC or LM in character. 

There is, however, a plausible explanation: that straight perturbation theory in U 
about the non- interacting limit is intrinsically inapplicable for | < r < 1. We believe 
this to be the case, for we have shown in section 3 that for [/ = the ground state 
is SC for all r < 1, and LM for r > 1. By contrast, finite-C/ NRG studies [12,13] 
yield a LM ground state for r > |, with no indication of a transition between LM and 
SC states at any finite interaction strength. For | < r < 1, therefore, the transition 
between SC and LM states occurs 'at' U = itself, and since the non-interacting and 
U > ground states are of different symmetry, one anticipates a breakdown of naive 
perturbation theory about the non-interacting limit. For r < | or r > 1 by contrast, 
the U = ground states are the same as those for U > (at least for sufficiently small 
U in the case of < r < |). Hence, as found in sections 5.1 and 5.2, one expects low 
order perturbation theory in U to be applicable over some finite-C/ interval; although 
we naturally anticipate such to become 'dangerous' asr— orr— > H-. 

6. Conclusion 

A summary of the present work is aptly illustrated by figure 7 which, for the symmetric 

soft-gap Anderson model, shows schematically the phase boundaries between strong 

1 

coupling (SC) and local moment (LM) phases in the U = U / Aq"*^ vs. r plane. In the 
non- interacting limit, we have shown in section 3 that for r < 1 the ground state is SC, 
while for r > 1 it is LM. From finite-C/ NRG studies [12,13] by contrast, it is known 
that the ground state is exclusively LM for r > |. In consequence, for r > 1, second 
order perturbation theory in U about the non-interacting limit as considered in section 
5 appears able to capture at least the intial evolution of the LM state upon increasing 
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U from zero, yielding in particular the characteristic low-frequency spectral signature of 
the LM phase found in NRG studies [12], viz D(uj) ~ |cj|''. For ^ < r < 1 by contrast, the 
[/ = and U > ground states are fundamentally distinct, and we believe a finite-order 
perturbative treatment about the non-interacting limit to be simply inapplicable. 

Perhaps the most subtle regime is r < ^, encompassing as a particular case the 
'normal' Anderson model, r = 0. In section 4, general conditions were established for a 
SC state to be realized for [/ > 0. Physically, these amount to interactions having no 
influence in renormalizing the lowest-frequency asymptotic behaviour of the impurity 
spectrum D{uj), whose behaviour as a; —> is that of the non- interacting limit, viz 
Diuj) ~ \oj\~'^ — producing the characteristic spectral signature of the SC state observed 
in NRG studies [12]. Significantly, one deduces in consequence a pinning condition 
(equation (4.7)) upon A{uo) = \ui\^D{uj), whereby A{uj = 0) at the Fermi level is pinned 
at its non-interacting value for all U and r where a SC state obtains, and which represents 
a natural generalization of the familiar pinning condition D{uj = 0) = I/ttAq V U 
characteristic of the r — Anderson model. At the level of second order perturbation 
theory in U (section 5.2), such a state of affairs is realized in practice for r < |, where the 
characteristic spectra A{u!) (figure 5) are strikingly reminiscent of their counterparts for 
r — (figure 6), with evolving Hubbard satellites upon increasing U and the emergence 
of a low-frequency Kondo scale refiected in the progressively narrowing width of the 
generalized Abrikosov-Suhl resonance in A{u!). 



U 




Figure 7. Schematic phase boundaries between SC and LM states in the U = U/Aq 
versus r plane; details in text. 



Figure 7 also shows clearly the limitations of finite order perturbation theory in 
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the interaction strength. For r < |, NRG studies yield both SC and LM phases 
[12,13], the critical line Uc{'^) diverging as r — >^ 0, reflecting the fact that the r = 
Anderson model is a Fermi liquid for all finite U; and vanishing as r — | (although 
whether it does so continuously as r |— we regard as not wholly settled). It is 
in part this feature — the existence of a transition between a LM state and a SC (or 
generalized Fermi liquid) state — that renders the problem of generic interest. But low 
order perturbation theory, capable though it is of describing the initial evolution of the 
SC state upon increasing U from zero, will naturally produce a SC state for all U and 
cannot therefore delineate its boundary. To deal with this central issue, as well as to 
describe successfully the regime ^ < r < 1, one requires an inherently non-perturbative 
theory that is capable of capturing both LM and generalized Fermi liquid phases, and 
hence the transition between them. For this, we believe a rather radical departure from 
conventional theoretical approaches is required, and will turn to one such in a subsequent 
paper [14]. 
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